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In this paper we examine the extent to which black hole quasi-normal modes (QNMs) could be 
used to test the no-hair theorem with future ground- and space-based gravitational-wave detectors. 
We model departures from general relativity (GR) by introducing extra parameters which change the 
mode frequencies or decay times from their values in GR. With the aid of Bayesian model selection, 
we assess the extent to which the presence of such a parameter could be inferred, and its value 
estimated. We find that it is harder to measure the departure of the mode decay times from their 
GR values than it is with the mode frequencies. The Einstein Telescope (ET, a third generation 
ground-based detector) could detect departures of as little as 8% in the frequency of the dominant 
QNM mode of a 500 black hole, out to a maximum range of ~ 6Gpc {z ~ 0.91). The New 
Gravitational Observatory (NGO, an ESA space mission to detect gravitational waves) can detect 
departures of ~ 0.6% in a 10* Mq black hole to a luminosity distance of 50Gpc {z ~ 5.1), and 
departures of ~ 10% in a 10^ Mq black hole to a luminosity distance of ~ 6 Gpc. In this exploratory 
work we have made a specific choice of source position (overhead), orientation (inclination angle of 
7r/3) and mass ratio of progenitor binary (mi/m2 = 2). A more exhaustive Monte Carlo simulation 
that incorporates progenitor black hole spins and a hierarchical model for the growth of massive 
black holes is needed to evaluate a more realistic picture of the possibility of ET and NGO to carry 
out such tests. 

PACS numbers: 



I. INTRODUCTION 

Merging compact binaries consisting of neutron stars 
or black holes are promising sources jlj of gravitational 
radiation for advanced gravitational wave detectors that 
are currently being built [2 [3] . These sources are also 
potential laboratories for testing general relativity (GR) 
in the strong field regime [IHH]- The result of such an 
event (with sufficiently massive progenitors) is a highly 
perturbed black hole (BH) which rapidly returns to its 
quiescent state (i.e., a Kerr BH [9 ), through the emis- 
sion of gravitational waves (GWs). The perturbation, 
and therefore the emitted radiation, can be expanded in 
the natural basis of oscillations, the quasi-normal modes 
(QNMs) [TU] (for a recent review on QNM see, e.g., [TT]). 
A Kerr black hole is characterised only by its mass and 
angular momentum, and so are the complex frequencies 
of its QNM oscillations [TU [13j, although the relative 
amplitudes of the modes depend on the specific details 
of the excitation. 

Detection of the characteristic ringdown GW signal of 
a BH would, therefore, allow a direct test of the no-hair 
theorem [Hj, and hence GR, through the comparison of 
frequencies and decay times of these modes with the pre- 
dictions of GR for a BH with certain mass and spin. 
In practice, the detection and discrimination of multiple 
modes is essential, as it is first necessary to infer the mass 
and spin of the black hole before checking for consistency 
between the modes. If any of the modes have some pa- 
rameter dependence, other than mass and spin, then the 
mass and spin obtained from these modes will not be 
consistent with that obtained from the others, and thus 
the source of emission must not be a Kerr black hole. 

In Ref. [13] Dreyer et al first developed the formalism 



for testing GR with QNM. They also suggested the test 
of the no-hair theorem through the measurement of more 
than one mode. Berti et al, [TSl [iQ investigated the ac- 
curacy of measurement of individual mode parameters 
using a Fisher matrix analysis and estimated the resolv- 
ability of individual modes in the complete signal as a 
function of signal-to-noise ratio. They conclude that the 
presence of a second mode can be inferred as long as 
the signal-to-noise ratio (SNR) is larger than a critical 
value, under the assumption that the presence of a ring- 
down signal has been confirmed and the parameters of 
the dominant modes are reliably measured. The critical 
SNR depends on the mass ratio of the progenitor binary 
but an SNR of 20 should suffice if the mass ratio of the 
progenitor binary is q — mi/m2 2. Kamaretsos et al 
showed recently that using black hole ringdown signals 
that are emitted after a binary merges, it might be pos- 
sible to recover the mass ratio of the progenitor binary 
from relative amplitudes of the QNMs fTT] . Their signal 
model, however, is based on non-spinning black hole bi- 
naries. Spins might induce systematic effects that could 
make it difficult to reliably measure the mass ratio of 
progenitor binary. However, we have to await ringdown 
models for spinning black hole binaries to assess what 
might or might not be possible to measure. 

Kamaretsos et al [T7] also proposed different ways of 
testing GR using QNMs. In particular, they proposed 
that it is sufficient to use any three of the observed fre- 
quencies and decay times for checking the consistency of 
the ringdown signal with GR. This minimal set of param- 
eters could make the test far more effective than trying to 
resolve the different mode frequencies and decay times. 
This is one approach we shall follow in this work. 

More specifically, in this paper we explore the fea- 
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sibility of the proposal in Ref. [T7] for a limited set of 
sources and parameter values. Our study is also limited 
to sources with mass ratio 1:2. We find that for the spe- 
cific ringdown model that neglects progenitor spins, it is 
possible to infer the total mass and spin of the black hole 
and the mass ratio of the progenitor binary from the ob- 
served ringdown GW, using a phenomenological mapping 
from these to the mode parameters. Having established 
this, we extend our waveform model to include arbitrary 
deviations from the predicted values of the frequencies 
and decay times, which can be estimated either along 
with the physical parameters of the source or on their 
own, if we assume no relationship between the complex 
QNM frequencies and the mass and angular momentum 
of the BH. This gives us two methods of testing GR: 

1. Measure the QNM parameters individually and 
check for consistency between the implied mass and 
angular momentum from each of the modes. 

2. Assume that the GR waveform is broadly correct 
but allow deviations in the parameters of a subset 
of the QNMs away from their predicted values and 
then perform model selection between the GR and 
non-GR models. 

We study both of these techniques in the context of fu- 
ture space- and ground-based detectors, particularly the 
Einstein Telescope [181 US] and a rescoped version of the 
Laser Interferometer Space Antenna called New Gravita- 
tional Observatory (NGO). 

The rest of the paper is organised as follows: In Sec- 
tion HJ we define the waveform model and describe the 
analysis methods used in testing the no-hair theorem. In 
Section |III| we will discuss the expected sensitivities of 
ET and NGO and the distribution of the signal-to-noise 
ratio of ringdown signals in these detectors. Section |IV| 
describes the simulations of the tests of the no-hair the- 
orem using two different methods: (1) consistency of the 
various parameters characterizing a ringdown signal and 
(2) Bayesian model selection applied to simulated ring- 
down signals buried in Gaussian background. Section [V| 
contains comparisons of the two methods used in testing 
GR and our conclusions. 



II. SIGNAL MODEL AND ANALYSIS METHOD 

In this Section we will discuss the nature of the ring- 
down waveform used in this study. The main focus will 
be to use a superposition of quasi-normal modes con- 
taining not only the dominant mode but also the first 
two sub-dominant ones. Ringdown signals that we study 
are assumed to be emitted by deformed black holes that 
form from the merger of compact binaries consisting of 
non-spinning components in quasi-circular orbits. The 
latter assumption allows us to use a phenomenological 
waveform model based on numerical relativity simula- 
tions; the observed signal is characterised by only three 



intrinsic parameters: the mass and spin of the final black 
hole and the mass ratio of the progenitor binary. 

A. Ringdown model 

A perturbed black hole emits a spectrum of modes 
characterized by three numbers {I, m, n). Indices I = 
2,3,..., and m = —I, . . . , +1, are the well-known spher- 
ical harmonic indices, and n — 0, 1, 2, ... , is the mode 
overtone index. Overtones other than the fundamental 
n — mode are not excited with significant amplitudes 
and have much shorter damping times |15) . We shall 
therefore only consider the fundamental mode and drop 
the index from further discussion. 

The two polarisations of the gravitational waveform, 
/i+ and hx , emitted by a BH of mass M during its ring- 
down are described as the sum over the QNMs, 

l,m>0 

(1) 

A'l 

hx{t) = -— ^;|m|e-*/"'"r>^"sinK„t-m0), 

l,m>0 

(2) 

for i > and h+ = = 0, for t < 0, t ~ being 
the start of the ringdown signal. Here r is the lumi- 
nosity distance to the black hole, Ai\^\(q) are the mode 
amplitudes that depend only on the ratio q — toi/to2 
(toi > 7712) of the component masses of the progenitor 
binary, and Tim(M,j) and u!im{M,j) are the character- 
istic mode damping times and frequencies. The mode 
damping times and frequencies depend only the black 
hole mass M and its spin j, but they are difficult to 
compute analytically. It is necessary to use numerical 
methods to compute them jTT] . Fits to some of the lower 
order modes can be found in Berti et al [15]. t S [0, tt) is 
the angle between the BH spin-axis and the line-of-sight 
to the observer and G [0, 27r] is the azimuth angle of 
the black hole with respect to the observer. F[.™ (t) are 
the sum of spin —2 weighted spherical harmonics |16j . 

yI'^(l) = _2l"'"(i, 0) + (-1)' -2Y'-"'{l, 0), (3) 
^^"(0 = -2i"'"(t, 0) - (-1)' -2Y'-"'{i, 0). (4) 

Although the waveform in Eqs. ([I]) and ^ contains 
a summation of modes over all values of I and to, the 
relative amplitudes of the higher order modes are signifi- 
cantly less than those of the lower, so we restrict the sum 
to the most significant modes. The waveform considered 
for this analysis is a superposition of the rt = 0, / = 2, 
m — 1 and n = 0, Z = m = {2,3,4} modes (which will 
hereafter be referred to as the 21, 22, 33 and 44 modes, 
respectively) with the 22 mode being dominant. 

The relative dominance of the higher order modes is de- 
pendent on the specific excitation that occurred (through 
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the form: 
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FIG. 1: Fits to the amplitudes of modes excited in the process 
of a binary black hole merger as a function of the symmetric 
mass ratio of the progenitor binary. 



Aim) and the position of the observer (through Yim)- 
This depends on both the mass ratio q = mi/ni2 > 1 of 
the compact binary components prior to merger, and the 
inchnation angle i. We used a phenomenological model 
based on fits to a set of numerical simulations to describe 
the mode amplitudes. We used the data in Kamaretsos et 
al |17) , but instead of using the mass ratio q we used the 
symmetric mass ratio J' = mi 7712/(7711+7712)^ = q/{l+q)'^ 
to derive new physically motivated fits. The 22 mode is 
expected to grow linearly with the symmetric mass ratio 
as noted in Refs. |5nj [21]. Modes with odd I or m are 
not expected to be excited when the progenitor binary 
is of equal mass and so we expect their behaviour close 
to = 1/4 to be some power of (1 — 4;^). These con- 
siderations led us to the following forms of the various 
amplitudes: 

^22(1^) = 0.864 z., (5) 
A2i{v) = 0.52(1- 4 7.)°-"i ^22(7^), (6) 
A^^{y) - 0.44(1-4 7.)°-^^ A^^i^). (7) 
^44(1.) = [5.4(7^-0.22)^0.04] ^22(7^). (8) 

Figure [T] plots the above fits together with the am- 
plitudes of the various modes derived from numerical 
simulations. The current data seems to indicate that 
A44 = A4^i{v) / A22{v) has a minimum at v = 0.22. How- 
ever, the amplitude of the 44 mode derived from numer- 
ical simulations is not very reliable and so this might be 
just an artifact of bad data. More accurate simulations 
are needed to validate this result. 

To complete the signal model we have to specify the 
damping constants T£m- Instead of damping times r^m 
it is customary to use quality factors Qim defined by 
Qim = ^tmJiml^- Bcrti et al |15j provide a simple fitting 
formula to many mode frequencies and quality factors of 



Q = qi + q2{l - jY' 



(9) 
(10) 



where /i, /2, /s, qi, (?2 and q^ are fitting constants. For 
the modes considered in this paper Table|T]lists the fitt mg 
constants. 



(/,m) 


/i 


/2 


h 




92 


<i3 


(2, 1) 


0.6000 


-0.2339 


0.4175 


-0.3000 


2.3561 


-0.2277 


(2, 2) 


1.5251 


-1.1568 


0.1292 


0.7000 


1.4187 


-0.4990 


(3, 3) 


1.8956 


-1.3043 


0.1818 


0.9000 


2.3430 


-0.4810 


(4, 4) 


2.3000 


-1.5056 


0.2244 


1.1929 


3.1191 


-0.4825 



TABLE I: The fitting constants in Eqns. Q and ([To|, for the 
21, 22, 33 and 44 modes, for the dominant n = overtone 

m. 



The gravitational-wave strain observed by a detector 
due to the excitation of QNM can be expressed as 



h = F+liA 



(11) 



where and -Fx are the detector antenna pattern func- 
tions, which depend on the position (0, (p) of the black 
hole relative to the detector and polarisation angle of 
the radiation (see Ref. [T7] for details). The measured 
strain, therefore, depends on the following 9 parameters: 



= {Af , I/, r, e, ip, tp, b, (j), to} , 



(12) 



which includes the epoch to when the signal arrives at 
our detector^. 

Figure[2j left panel, shows the response of a signal con- 
taining all four modes considered in this work as well as 
the phase evolution of each of modes 21, 22, 33 and 44 for 
a 5 X 10^ A/0 black hole at a distance of 1 Gpc, assumed 
to be formed from a binary of mass ratio q = 2. The 
black hole is assumed to be overhead with respect to the 
detector, the orbital plane making an angle of t = 7r/3 
radians — a sub-optimal orientation but one for which 
the sub-domimant modes will have a nonzero amplitude. 
The polarization and azimuth angles are randomly cho- 
sen to be 7/' = 2.67 radians and (j) = 2.31 radians. 

The overall amplitude is very nearly the same as that 
of the 22 mode. The amplitude of the 33 mode, which 
is the second most dominant after 22, is about a third 
of the 22 mode, followed by the 21 mode, which has a 
slightly smaller amplitude. The 44 mode, whose ampli- 
tude is about 12% of the 22 mode, has negligible effect on 
the overlall signal. Although 21 and 33 are of roughly the 
same amplitude, the energy in 33 mode, whose frequency 



^ Note that Eqs. ljT| and are written assuming to = 0. A 
nonzero value of to can be trivally included in the signal model 
by writing t as t — to ^^nd assuming /i+,x (i) = for t < to- 
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is roughly twice that of 21, will be significantly larger 
as we will see in Sec. Ill while discussing the relative 



signal-to-noise ratios of the different modes. The inter- 
ference between the different modes causes features that 
become apparent in the spectrum of the modes shown 
in the left hand panel of Figure [2j which we will discuss 
further in Sec. IIIII 



T-L, p{0\'H) is the prior distribution of the parameters 
given the signal model and p{d\9,T-L) is the likelihood of 
the data for a particular set of parameters 9 : 



\d^-Hff,0)\^^ 



(17) 



B. A generalised model for the QNM 

In order to test GR, we extended the waveform model 
so that the frequencies and decay times of the modes 
were allowed to be dependent not only on M and j, but 
also other dimensionless parameters. More specifically, 
we considered that frequencies ujim depended on three 
parameters (M, j, Acdim) and damping times r/„i also de- 
pended on three parameters (M, j, Af/m). Furthermore, 
we assumed that the dimensionless parameters Aujim and 
Ar/m were independent for each mode. Following Li et al 
, in this generalised model, the frequencies i^;m,nonGR 
and decay times T;m^nonGR were expressed as 



Wfm^nonGR 
'Hm.nonGR 



l^im,GR (1 
T;m,GR (1 - 



Af;„) 



(13) 
(14) 



where uJim,GR and Thn,GR are the frequencies and decay 
times of modes as in GR. The signal produced by the GR 
hypothesis is a special case of the generalised model, in 
which AcD;„i = Afjm = for all I, m. 



C. Bayesian analysis 

Having described the waveform model and its parame- 
ters (contained in a parameter vector 9), we will now de- 
scribe how these parameters are estimated from data con- 
taining a ringdown signal. We assume that the data from 
the gravitational wave detector in the frequency domain d 
contains both the ringdown signal h{ f; 9 ) and some addi- 
tive Gaussian noise with known power spectrum Sh{f)- 
Thus, the data d is assumed to be di — h{fi;9) + fii, 
where i is the index of the frequency bin. The noise 
power spectra used for ET and NGO are given in Sec- 
tion [III A[ As we perform our analysis in the frequency 
domain, we use the Fourier transformed signal model 
h[f) = h{t)e-^'''f^dt computed with the FFTW 
package. 

Our goal is to compute the posterior probability dis- 
tribution (PDF) (see, for instance, Ref. 123; ), of the pa- 
rameters ^(6*1^, v.), 



p{9\d,n) 



p{d\9,H)p{9\H) 
p{d\H) 



(15) 



where p{d\'H) is the evidence, or marginal likelihood, of 
the model 

p{d\H)^ I p{9\H)p{d\9,H)d9, (16) 



Posterior distributions for particular parameters of inter- 
est, e.g. the AtD/m and Afim, are computed by marginal- 
ising the PDF over all other parameters. We also com- 
pute the Bayes factor Bi^j between various hypotheses, 
Bi^j = p{d\i) / p{d\i) , which is the evidence ratio between 
hypotheses i and j. Due to the large range of this quan- 
tity we will always use the natural logarithm of this, 
logSij , in section 



IV C 



To compute the PDFs and evidence we use the 
inspnest implementation of the nested sampling algo- 
rithm, described in j23j . modified to use the ringdown 
signal model, which is available as part of the LSC Al- 
gorithm Library [24]. The end product of the analysis is 
the model evidence and a set of samples from the poste- 
rior PDF that are histogrammed to produce the figures 
below. 



III. VISIBILITY OF THE QNM SIGNALS 

In this Section we will explore the visibility of quasi 
normal mode signals. We will begin with the expected 
sensitivities of ET and NGO and then go on to explore 
the visibility of the signal in these two detectors. The 
signal-to-noise ratio (SNR) of a ringdown signal depends 
quite critically on a number of intrinsic and extrinsic pa- 
rameters of the source. We will vary most of these pa- 
rameters to obtain a distribution of the SNR. 



A. Sensitivity curves 

In this paper, the ET and NGO detectors are consid- 
ered. For simulations concerning ET, the power spec- 
tral density corresponding to the ET-B sensitivity curve 
is considered, described by Sh{f) ~ 10~^"/i„(/)^ Hz^^, 
where h„{f) is given by 

h,,{f) = 2.39 X 10-27a;-i5-«4 ^ 0.349a;-2 "5 



1.76x" 



0.409a;' 



and X — //lOOHz. For NGO, the sensitivity curve asso- 
ciated with the L — Ix 10^ m arm, 4-link mission studied 
in Ref. [25 , which corresponds to a power spectral den- 
sity given by 



Shif) 




2Lf 



0.41c 



i4.Sacc + So) , (18) 
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FIG. 2: Left: Strain amplitude of a quasi-normal mode signal from a black hole that forms from the merger of a binary of 
(observed) total mass 5 x 10^ Mq and mass ratio q — 2 sd, 1 Gpc. We have plotted the first four dominant modes, 21, 22, 33 and 
44 together with their superposition. Right: The signal-to-noise ratio integrand of the same signal dp'^/df = \H{f)\^ /Sh{f), 
where Sh{f) is taken to be that of NGO. The presence of the 33 and 44 mode can be clearly seen in the overall spectrum, 21, 
however, is buried under 22. 




FIG. 3; Distribution of the signal-to-noise ratios in different quasi-normal modes (21, 22, 33 and 44) for sources located at 
random positions on the sky, with random inclination and polarisation angles. The left plot is for QNM resutling from the 
merger of a 500 solar mass binary observed in ET and the right plot is for a 5 million solar mass binary observed in NGO. In 
both cases the mass ratio of the binary is assumed to be g = 2 and the source is assumed to be at 1 Gpc. 



where. So = 1.153 x lO"^^ j^2jj2-i and 



Sarr. = 1-37 X 10 



-32 



lO-^HzX /27r/ 



IHz 



(19) 



B. Visibility of ringdown signals 



power of the test. The SNR depends both on the in- 
trinsic parameters of the black hole, its mass and spin 
angular momentum, as well as its extrinsic parameters, 
its distance and various angles describing the orientation 
of its spin axis, its location on the sky and the polarisa- 
tion of the radiation. In this Section we will look at the 
distribution of the SNR in ET and eLISA/NGO. 

The SNR of a ringdown signal observed using matched 
filtering is given by 



An important requirement for tests of the no-hair theo- 
rem is good visibility of the signal in our detector: larger 
the signal-to- noise ratio (SNR), greater would be the 
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d/, 
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FIG. 4: Projections in the (M, j)-plane of the 90% confidence hmits on CJ22, T22 and tJss (blue, blue dotted and red lines 
respectively) for injections of signals consistent with GR for M = 500 M© (top-left at 125 Mpc; SNR = 2 888), M = 1000 M©, 
(top-right at 225 Mpc; SNR = 2 423); M = lO'^ Mq (middle-left at 125 Mpc; SNR = 63), M = 10^ A/q (middle-right at 125 Mpc; 
SNR= 1756), M = 5x IO'^Mq (bottom-left at IGpc; SNR = 6 377) and M = 10® Mq (bottom-right at 1 Gpc; SNR = 115154). 
The injected value is denoted in each case by a diamond. 
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where Sh{f) is the detector noise power spectral den- 
sity and H(f) is the Fourier transform of the signal in 
Eq. ( |11[ ). The limits on the integrand can, in principle, 
be from to oo. However, in order to prevent the "junk" 
radiation that occurs in the Fourier transform due to 
abrupt cutoff of the time-domain signal from corrupt- 
ing the SNR we choose finite values for both them. More 
concretely, we choose /low to be half of the smallest sig- 
nal frequency (which in our case is that corresponding 
to the 21 mode) and /high to be twice the highest signal 
frequency (which in our case is that corresponding to the 
44 mode). 

p is the SNR of the signal containing all the modes 
considered in this paper. To test the no-hair theorem 
it is essential that the SNRs of the sub-dominant modes 
(21, 33, 44) are comparable to the dominant 22 mode. To 
assess their importance, we shall also separately consider 
the SNR pfm of each mode: 



2 
Pin 



he 



Shif) 



df, 



(21) 



where Hf, 



is the Fourier transform of the signal in 
Eq. (Ill but with and hx containingly only the rel- 



evant mode. We will compoare the distribution of pim 
for £m = 21, 22, 33 and 44 modes but it should be kept 
in mind that p^ is not the sum over different p^^. Even 
so, the relative strengths of different modes gives us an 
indication of how good we can expect our tests of the 
no-hair theorem to be. 

Figure [2j right panel, plots the SNR integrand dp^ /df 
of the various modes. The overall signal gets most of its 
contribution from the 22 mode but other modes signifi- 
cantly alter the phasing of the signal (left panel) and its 
spectrum. The overall signal has bumps caused by the 33 
and 44 modes but the 21 mode, whose frequency is close 
to that of 22, has little effect on the overall amplitude or 
the phasing of the waves. Note that the spectrum of 33 is 
far larger than the 21 mode although they are of similar 
amplitude. 

Figure [2] gives a qualitative understanding that the 
subdominant modes have significant effect on the phas- 
ing and spectrum of the signal. But we need to evaluate 
the distribution of the SNR to better gauge the relative 
importance of the different modes. The response of our 
detectors to a ringdown signal depends on a total of nine 
parameters as enumerated in Eq. (12). To compute the 



distribution of the SNR, however, we don't need to vary 
all nine parameters. We know that the SNR depends 
inversely on the luminosity distance r of the source and 
so this can be fixed to any value we wish and we take 
it to be 7' = 1 Gpc. The distribution does not depend on 
the epoch of the signal tg nor on the phase (f). We should 
study the distribution as a function of M and v while 
varying the remaining four parameters. However, since 
our goal is to explore the test of the no-hair theorem for 
a small sample of signals with high SNR we will restrict 
ourselves to only a couple of sources. We will consider a 
BH that results from the merger of a 500 Mq binary in 



the case of ET and a 5 x 10^ M© binary in the case of 
NGO. In both cases we take the symmetric mass ratio to 
be g = 2 or 1/ = 2/9. 

The cumulative distribution of the SNR is shown in 
Figure [3| We see that there is a 25% chance that the SNR 
will be greater than 300 in ET and 3000 in NGO. The 
corresponding SNR in the 33 mode is 90 and 900 for ET 
and NGO, respectively. Even the least dominant 44 mode 
has SNRs of 20 and 200. The SNRs drop off as inverse 
distance and at redshifts z ~ 1-3, from within which we 
can expect the event rate of massive black holes to be 
in excess of several per year, the SNRs will be still large 
enough that NGO should be able to carry out meaningful 
tests of GR for a large fraction of sources detected. Only 
very rare closeby events will allow such tests in the case 
of ET. 



IV. TESTS OF THE NO-HAIR THEOREM 

In this Section we will explore the ability of ET and 
NGO to test the no-hair theorem if they were to detect a 
ringdown signal from a black hole that resulted from the 
coalescence of a binary. We will begin with a summary 
of the parameters used in our simulations followed by 
a description of the two different approaches that were 
used to test the no-hair theorem. The first method uses 
consistency of the various mode frequencies and damping 
constants with GR predictions and the second method 
is a Bayesian model selection approach that addresses 
which of a class of different models best describes the 
underlying signal. We find that the latter approach is far 
more powerful in testing the no-hair theorem. 

The model selection approach followed in this paper 
is computationally pretty expensive. We have therefore 
selected a small sample of signals (two in ET and four in 
NGO) to assess how well future detectors are able to test 
the no-hair theorem. The SNR study from the foregoing 
Section can be used to conclude how effective such tests 
are for a random signal detected in our instruments. 



A. Choice of injection parameters 

The GW signal emitted from the QNM of a black hole 
ringdown as observed in a detector depends on the mass 
and spin of the final BH, the mass ratio of the progenitor 
binary and other extrinsic parameters that describe the 
orientation of the black hole, its distance from the detec- 
tor and its position on the sky. For the different sensitiv- 
ity bands of ET and NGO we chose the following range of 
source parameters. For ET, black holes of observed mass 
500 Mq and 1000 M© at luminosity distances, I?^, from 
125 Mpc and 225 Mpc respectively, out to 6.63 Gpc (cor- 
responding to redshift z ~ 1 [26]), were considered. For 
NGO, black holes of observed mass 5 x 10*^ and 10^ M© 
at Di, = 1-59 Gpc (with the upper limit corresponding 
to a redshift z ~ &), and black holes of observed mass 10^ 
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FIG. 5: Projections in the (M, j)-plane of the 90% confidence limits on a;22, T22 and 0^33 (blue, blue dotted and red lines, 
respectively) for non-GR injections of M = 500 Mq (top at 125 Mpc; with Alj22 = -0.01, SNR = 2 867 (left) and Aa>22 = -0.05, 
SNR = 2779 (right)), M = 10^ (middle at 125 Mpc; with Aa)22 = -0.01, SNR = 1753 (left) and Ati22 = -0.05, SNR = 
1735 (right)) and M = 10* (bottom at 1 Gpc; with Aci22 = -0.001, SNR = 115130 (left) and Ad)22 = -0.005, SNR = 
115 031 (right)). The injected value is denoted in each case by a diamond. 
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and 10^ M© at Dl = 125 Mpc - 6.63 Gpc, were consid- 
ered. In both cases, the source position (a, (5) was set to 
be directly above the detector at the time of observation. 

For all systems, the inclination and polarisation an- 
gles (t,?/") were set to be (f ,f) and the azimuth was 
taken to be zero. For the first four systems described, 
the mass ratio q of the binary system prior to merger 
and the spin j of the black hole after merger were set 
to g = 2 and j = 0.61, respectively. For the final two 
systems outlined, however, q = \Q and j = 0.26, to rep- 
resent a more typical NGO binary system. For the sake 
of using the phenomenological fit in table |T] the binary 
components were assumed to be non-spinning. It is im- 
portant to note, especially for NGO, that the observed 
mass of the system is redshifted to be greater than the 
intrinsic mass of the system by a factor of (1-1- z). Here we 
report this observed, redshifted mass for injections and 
recovery. 

For injections consistent with GR, parameters AcDim = 
— 0, whereas for non-GR signals, either Aw22 or 
/S.T22 was varied in the range —0.01... — 0.1, with all 
other Awjm = ^hm = 0. 



B. Constraining QNM parameters 

Having established the parameters of our test sources, 
we now move on to consider the two tests of GR. The 
first method broadly follows the outline of Ref. [17 to 
estimate the parameters of each QNM in the ringdown 
signal. By choosing two of the measurements one can in- 
fer the true mass and spin of the black hole by excluding 
regions in the M, j plane which are inconsistent with the 
measured parameters. The third measurement is then 
used to confirm the consistency of the inference. In the 
case of GR the allowed area will intersect both the pre- 
vious areas as in Figure |4] Alternatively, if the signal is 
inconsistent with GR the intersections of the confidence 
regions will not agree, as in Figure [5] 

In order to compute the constraints, wc find the upper 
and lower 90% probability interval of the marginal poste- 
rior probability distribution of the parameters Ad);™ and 
Afim- In each case the injected values of these param- 
eters were set to 0, corresponding to the GR waveform. 
The priors for Af/m and AcDj^ were all uniform in the 
interval [—1.0, 0.3] for all analyses, and the M, j and q 
parameters were set to their injected values. As we are 
estimating the mode parameters directly, through Acjjm 
and Afjm , the values of M and j parameters do not play 
any role in the results other than to set the centre of the 
prior range. The reason for the asymmetric choice of in- 
terval is that by allowing Aa)22 to be greater than 0.3, 
the frequencies of the 22 and 33 modes collide, leading 
to strong correlations between the parameters of these 
modes upon recovery and confusion between the modes. 
Note that we allow the frequency of the 33 mode to vary 
downward into the band of the 22 mode, but this does 
not result in difficulties as the 22 mode is already within 



this band, making its identification easier. In all searches, 
the luminosity distance of the source and orientation an- 
gles are assumed known and fixed to the aforementioned 
values. 

We chose the Acj22, Af22 and ACJ33 parameters to 
perform our consistency test as these are the three re- 
covered parameters with the greatest precision for the 
signals we considered, and will therefore give the most 
stringent test. 

We first injected signals using the GR waveform and 
performed parameter estimation on the three test pa- 
rameters. The 90% probability limits on ujim and Tim 
were projected in the (M,j)-plane to show visually the 
consistency test between the three modes, and that they 
agree at the injected value (M, j) within the measure- 
ment accuracy. Figure [4] shows that for each system, the 
projections of w/m and r/,„ coincide at the same position, 
and the region of coincidence encloses the injected value 
of the mass and spin of the system, as expected. 

To contrast, in Figure [5] we show the corresponding 
plots where the injection is performed with deviations 
from GR of A(2;22 of -0.01 and -0.05 for both ET and 
the lO^M© NGO system, whilst for the IO^Mq NGO 
system, deviations of —0.001 and —0.005 were considered 
(all other parameters are taken to be the same). 

This demonstrates the feasibility of the method for 
testing the no-hair theorem using consistency of the 
modes, when the signal is strong and the modes are 
clearly distinguishable. Of course, we may not be fortu- 
nate enough to observe BH ringdowns from such nearby 
sources, in which case we expect the power of the method 
to diminish. 

To investigate the accuracy to which the mode param- 
eters are resolvable we performed a set of injections at 
luminosity distances spanning the entire distance range 
quoted in Section |IV A[ and estimated the AJ>22, A(2;33 
and Af22 parameters. For each black hole system con- 
sidered, the width of the 90% confidence intervals for the 
extracted values of Ail;;,„ and Af/™ were plotted against 
luminosity distance in Figure |6j for injections with GR 
waveforms. Our results are in agreement with similar 
studies carried out in Ref. [27] . From Fig.[6j it can be 
seen that the width of the 90% confidence intervals for 
the values of Aa;22,33 and Af22 increases with distance, 
as expected. Aa}33 is extracted with considerably less ac- 
curacy than Ac2;22, because for a black hole system with 
q = 2, as considered here, the 33 mode is significantly 
less excited than the dominant 22 mode, and thus the 33 
mode has a much lower SNR, resulting in poorer resolu- 
tion. By the same token, Af22 is extracted with consid- 
erably less accuracy than both Aw22,33, as there is less 
extractable information from the mode decay times as op- 
posed to the mode frequencies. In general, however, the 
relative weights and SNRs of each mode depends on the 
source parameters including orientation, and the noise 
curve of the detector. 

We also see that ET can resolve AJ)22 with errors of 
less than 1% at Dl = 4 Gpc for the systems considered. 
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FIG. 6: Width of the 90% confidence intervals for Acj22,Aii33 and Af22 (blue, red and blue dotted lines respectively) against 
luminosity distance for injections of 500 (top-left), 1000 (top-right), 10* (middle-left), 10® (middle-right), 5 x 10® (bottom- left) 
and IO^Mq (bottom-right) for ET (top) and NGO (middle and bottom) simulations. The limiting values in the IO^Mq plot 
are due to the posteriors expanding to fill the prior distributions when they are no longer resolvable. 
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in addition to resolving Aw22,33 and ^ti2 with errors 
of 1.5%, 2.5% and 10% for a black hole of 500 Af© at 
Dl = 6.63 Gpc. 

NGO can resolve ACD22 to ^ 0.6% accuracy at Z?^ = 
59 Gpc (z = 6) for both 5 x lO^Af© and 10^ sys- 
tems, showing good results across the range of masses de- 
tectable. Af22 (Acjas) can be resolved to ~ 3.5% (1.5%) 
and - 1% (0.15%) accuracy for 5 x 10^ Mq and 10*^ Mq 
systems respectively, at this distance. This trend is mir- 
rored for the 10® Mq system, with NGO able to resolve 
Ad)22 to ~ 2% accuracy at = 6.63 Gpc (z = 1), and 
A'r22 (AW33) to 10.5% (5%) at this distance. For the 
10^ Mq system however, resolution of Aa}22 is poor even 
at distances of Dl = 1 Gpc, suggesting that it would not 
be possible to conduct these discriminatory tests using 
such systems, as they are not expected to exist at dis- 
tances closer than 1 Gpc. 

The method outlined above to test the accuracy with 
which we can extract the three QNM parameters was 
applied to the case of non-GR signals (Acj/^ or Af;^ 
set to some value within the range quoted in Section 
IV A with all other Aw/^ = Afjm — 0) from systems 
of M = 500 Mq and 10^ Mq in ET and NGO. The QNM 
corresponding to these systems lie closest to the most sen- 
sitive frequency regions of the corresponding detectors, so 
this test is an optimistic case. Also, the distances used 
were from the lower end of the range quoted in Section 
IV A For each non-GR injection, the 90% probability 
limits on the values of Ac2;22, Acjsa and A'r22 were ex- 
tracted and used to project the corresponding confidence 
limits on uJim and t;,„ in the (Af, j)-plane. 

Figure [5] shows that even for deviations from GR of 
as small as 1% (corresponding to AW22 = —0.01) for 
ET and the 10® Mq NGO system, and as small as 0.1% 
(corresponding to Aw22 = -0.001) for the 10^ Mq NGO 
system, the projection of W22 does not intersect consis- 
tently with the projections of T22 and ^33, which were 
unchanged from their GR predicted values. When an 
additional parameter dependence is introduced for one 
mode parameter, the mass and spin derived from this 
modified parameter will not be consistent with the mass 
and spin derived from the other unchanged parameters. 
For a deviation of 5% for ET and the 10® Mq NGO sys- 
tem, and 0.5% for the 10* Mq NGO system, the projec- 
tion of ^22 does not even touch the projections of T22 and 
W33, let alone intersect at the correct value. This shows 
that in extracting the individual mode parameters from 
a QNM signal and projecting them in the (M, j')-plane, 
it is possible to determine if one, or many, of the mode 
parameters are not consistent with GR and have some 
parameter dependence other than on the mass and spin. 



C. Discriminating between models 

Although we have demonstrated the possibility of per- 
forming a consistency test between the QNM parameters 
by plotting them on the {M,j) plane, this method does 



not provide a quantitative measure of the consistency as 
such. In this section we seek to provide such a quantifi- 
cation by performing model selection on two competing 
models: the standard GR model for which we estimate 
M, j and q; and the extended model where we estimate 
Af , j and q but allow the individual mode parameters to 
vary away from their GR predicted values through the 
AuJim and Afim parameters. We shall see that even in 
the event where a source is too distant to allow the con- 
sistency method discussed in the previous section to be 
applied, it could still possible be to perform model selec- 
tion to distinguish between GR and non-GR models. 
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FIG. 7: The differences in evidence between the GR and non- 
GR models against Aii22 and AT22 (blue and blue dotted 
lines respectively). (Top-Left) For injections oi M — 500 at 
125 Mpc in ET, with SNR between 2 664 and 2 888. (Top- 
Right) For injections of 10** in NGO at 1 Gpc, with SNR 
between 114 909 and 115154. (Bottom) For injections of 
lO^M© in NGO at 125 Mpc, with SNR between 1 706 and 
1757. 

Bayesian model selection is a method of determining, 
for a given signal, whether the GR or non-GR model is 
more likely by comparing the evidences for each. Note 
that unlike in the previous case where we estimated the 
QNM parameters directly, here we are estimating the 
physical parameters of the source along with deviations 
from the QNM parameters, where all QNM parameters 
are allowed to vary together. The reason for this differ- 
ence is that in estimating the QNM parameters directly 
we would not have a corresponding GR model with any 
free parameters unless we also searched over Af, j and q. 
Because in the non-GR case we are measuring five param- 
eters, including the best-determined ones, the resulting 
joint probability distributions are strongly correlated, al- 
though the presence of the 21 mode and the assumption 
that Ar33 = help to break the degeneracy. 

Here we compare evidences for both GR and non-GR 
signals injected into simulated ET and NGO data as de- 
scribed above, but calculate the evidence for each model 
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using the nested sampling algorithm. The ratio of these 
evidences, known as the Bayes factor, quantifies the sup- 
port for one model over the other provided by the data. 
Full computation of the Bayesian evidence automatically 
takes into account the difference in dimensionality of the 
parameter space, which penalises the more complex non- 
GR model accordingly, as we shall see. 

Figure [7] displays how the Bayes factor between the 
GR and non-GR models changes as the deviation of the 
injected signal from GR increases. In the case of ET, we 
used the signal from a 500 Mq black hole at at distance of 
125 Mpc, giving a signal to noise ratio between 2 644 and 
2 888. The NGO simulation used both a 10® Mq system 
at a distance of 1 Gpc, with signal to noise ratio between 
114 909 and 115154, and a 10^ Mq system at a distance 
of 125 Mpc, with signal to noise ratio between 1 706 and 
1 757. For ET and the 10^ Mq NGO system, we varied 
both Aci)22 and Af22 in the range —0.1 to 0, whereas 
for the 10® Mq NGO system, we varied Ad)22 and Af22 
in the range —0.01 to 0, and in each case computed the 
evidence for both GR and non-GR models. As the de- 
viation from GR increases, the evidence for the non-GR 
model increases whilst the evidence for the GR model de- 
creases, causing the Bayes factor (the ratio of evidences 
for the two models in log space) to increasingly favour the 
non-GR model, as expected. This is because as the devi- 
ation from GR increases, the data fits the GR model with 
less accuracy, but the parameter freedom of the non-GR 
model permits a consistency with the data. If we assume 
a threshold of log BcR^non-GR = —10 (corresponding to 
prior odds of in favour of the non-GR model), below 
which a deviation is considered significant, then with the 
source parameters above we are able to distinguish devi- 
ations at the level of 1% with ET and NGO for the most 
sensitive Aw22 parameter. 

As with the mode resolution of testing GR, the power 
of the model selection test decreases with increasing 
source luminosity distance since the signal visibility de- 
creases. The model selection method was carried out for 
the sources above at distances spanning the entire range 
outlined in Section [IV A[ and the Bayes factor between 
the GR and non-GR minimal set plotted as a function 
of distance for signals with different deviations from GR. 
Figure |8] shows that for signals with larger deviations 
from GR, the maximal set non-GR model is favoured 
out to much greater distances. We see also that for GR 
injections, the GR model is indeed preferred, since the 
simpler model fits the data adequately. When the dis- 
tance is increased, the size of this effect decreases as the 
posterior distribution for the modified parameters grad- 
ually expands to fill the prior range. The power of the 
model selection method relies on the ability to exclude 
parts of the parameter space of the more complex model. 

For the purposes of comparison, if we again impose 
a threshold on the Bayes factor of logBcji^non-GR — 
— 10, we see that the non-GR model is favoured above the 
threshold at all distances for the 10® Mq system with at 
a deviation of only 0.6%, whereas deviations as small as 



0.1% can be detected out to a distance of ~ 10.7 Gpc. For 
the 10^ Mq NGO system, the non-GR model is favoured 
above the threshold up to ^ 1 Gpc for a 2% deviation 
and ^ 6 Gpc for a 10% deviation. For ET, we are able to 
cross the threshold up to ^ 1.61 Gpc for a 2% deviation 
and 6.63 Gpc for an 8% deviation. 



V. CONCLUSIONS 

In this paper we have investigated how well quasi- 
normal modes could be used to test general relativity. 
Our work is based on a quasi-normal mode signal model 
|17) which assumes that the progenitor binaries are non- 
spinning. We have specifically looked at the ability of 
future ground- and space-based detectors in testing the 
black hole no-hair theorem. More specifically, we have 
investigated how well future interferometric gravitational 
wave detectors can test deviations of quasi-normal mode 
frequencies and damping from their general relativistic 
values. 

As expected, the / = 2, m = 2 mode of the signal 
is the most clearly measurable for the sources we con- 
sidered, which can be used to infer the mass and spin 
of the black hole. Measurement of a third parameter, 
i.e. the frequency of the Z = 3, to = 3 mode, is then 
used to confirm consistency between the modes. If the 
true signal is consistent with GR then the measurement 
of the different parameters would be consistent with one 
another. We find that a 10% deviation in 0^22 is clearly 
discernible for a 500 Mq source at 1.25 Gpc in ET, as well 
as 10^ and 10® Mq sources at 1.25 Gpc and 10 Gpc, re- 
spectively, in NGO. Within the reach of NGO, the event 
rate could be pretty high and so no-hair tests are promis- 
ing in this case. Binary coalescences of intermediate mass 
black holes within 1.25 Gpc are highly unlikely and hence 
ET might not be able to carry out such tests, except in 
the case of rare close by events. 



In Section IV C therefore, we applied Bayesian model 
selection to obtain a more robust and quantitative mea- 
sure of the consistency of the data with GR vs a gener- 
alised theory where the mode parameters depended on 
an extra parameter other than the black hole mass and 
spin (i.e., 'hairy' black holes). Using this technique, we 
are able to measure deviations at the 10% level in the 
W22 parameter out to ~ 6 Gpc for a 500 Mq source with 
the Einstein Telescope. With NGO, we are able to mea- 
sure deviations at the 10% level at 6 Gpc with a 10^ Mq 
source, and the 0.6% level at z ~ 5.1 with a 10® Mq 
source. 

In the above analysis, we have assumed that the lo- 
cation and orientation of the GW source is known prior 
to the analysis — a significantly simplified analysis com- 
pared to the full problem of determining these parame- 
ters alongside the test of relativity. Another limitation 
of our study is that it is restricted a small number of sys- 
tems with very specific source location, orientation and 
mass ratio of the progenitor binary. Moreover, due to the 
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FIG. 8: The differences in evidence for tfie GR and non-GR models is plotted against the luminosity distance to the source, 
for observation of ringdown signals from a 500 Mq black hole in ET (top panels), 10* and 10^ M© black hole in NGO (middle 
and bottom panels repectively). The non-GR model is simulated by varying either Aa;22 (left panels) or Af22 (right panels) 
as compared to their GR values of 0. Aa)22 and Af22 are varied over the range to —0.1 in steps of —0.01 in the case of 
ET and the lO'' A/q NGO system, and over the range to -0.01 in steps of -0.001 in the case of the 10® M© NGO system, 
while keeping all other parameters as for GR. The different curves in each panel correspond to different values of the non-GR 
parameter, starting with for the top most curve and changing by either —0.01 (top and bottom panels) or —0.001 (middle 
panels). A difference in evidence of —10 or smaller is considered good enough to discriminate a non-GR model from GR model. 
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lack of accurate models, the study has used a very sim- 
ple model of the ringdown signal that neglects the effect 
of initially spinning black holes of the progenitor binary. 
The presence of spins could significantly alter the spec- 
trum of the modes excited and their relative amplitudes, 
which could impact our ability to infer the mass ratio of 
the progenitor binary and the no-hair tests studied in this 
paper. Nevertheless, the methods we have developed in 
this work are equally applicable to the more general case, 
which should be performed as part of a follow-up study. 
Such a study will be helpful in a more robust evaluation 
of the potential of future detectors in testing the no-hair 
theorem. 
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